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Abstract 

Lattice QCD with additional chiral-invariant four-fermion interactions is 
studied at nonzero temperature. Staggered Kogut-Susskind quarks are used. 
The four-fermion interactions are implemented by introducing bosonic auxil- 
iary fields. A mean field treatment of the auxiliary fields is used to calculate 
the model's asymptotic scale parameter and perturbative thermodynamics, 
including the one-loop gluonic contributions to the energy, entropy, and pres- 
sure. In this approach the calculations reduce to those of ordinary lattice 
QCD with massive quarks. Hence, the previous calculations of these quan- 
tities in lattice QCD using massless quarks are generalized to the massive 
case. 
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I. INTRODUCTION 



One of the remaining challenges within the standard model of particle physics is to 
understand the high temperature behavior of QCD with light quarks. This includes deter- 
mining the location and order of the chiral phase transition, associated with the restoration 
of the SU(A^/)l eg) SU(A^/)r flavor symmetry, and whether it obeys the expected dimensional 
reduction scenario In the high temperature plasma phase, the status of the U(1)a sym- 
metry, the equation of state, and the gluon screening masses, necessary for high temperature 
perturbation theory, still remain to be determined. 

Lattice gauge theory appears to be an excellent nonperturbative tool for studying these 
issues. In fact, many aspects of the deconfining phase transition and the high temperature 
plasma phase in the gluonic sector are now well understood from Monte Carlo studies of 
pure SU(3) gauge theory 0. Unfortunately, simulations of full QCD, including dynamical 
quarks, are much more computationally intensive and so must be performed on smaller 
lattices with relatively large lattice spacings and/or finite size effects. Ref. contains some 
recent reviews. 

In an effort to more closely approach the true physics of continuum QCD, many re- 
searchers have turned to using improved or "perfect" lattice actions, which reduce 
discretization errors and lattice artifacts. A somewhat different approach is the xQCD 
model 0-0. In this model, extra auxiliary fields are introduced, which upon integration 
yield chiral invariant four-fermion interactions of the form 2jv^/4 [(^V')^ ~ (V'75''"V')^]- Such 
interactions are irrelevant to the long distance physics and do not survive in the continuum 
limit. However, when using the hybrid Monte Carlo algorithm the auxiliary fields entering 
the fermion determinant allow one to avoid the usual zero mass singularity when inverting 
the Dirac operator and perform simulations directly in the chiral limit, (i. e. the bare quark 
mass can be set exactly to zero.) 

To be specific, the xQCD action for Nf fiavors of staggered quarks on an anisotropic 
lattice with spatial (temporal) lattice spacing a (ot-) is 

'5 = E /5-E^u(^) + /5rE^o, + E E«Vx"(a;)Q(x,2/)x"(y) 



X 



a=l x,y 



+ ^7E«V[a^(5)+7r^(5)], (1) 



withQ 



3 



Q{x,y) = Y.-^j(^'y) +'yFMo{x,y) + 6^^y^ a[a{x) +ie{x)TT{x)] (2) 

j=l <x,x> 



^ The model as considered here has a U(1)l U(1)r global symmetry, generated by 1 and 75X3 
(t3 being the flavor equivalent of 75). One can also consider more physically realistic SU(A'^^)l <8> 
SU(A^/)r symmetries. However, for the calculations considered in this paper, only the total number 
of quark flavors is relevant. 
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and 7 = l/G. The staggered hopping term is 



Mu{x,y) = ^r]^{x)[U,,{x)6y^^+i) - U^{x - u)Sy^^_o , (3) 

where r]y{x) = (^—i"j^o+---+x^-i ^ rjj^g auxihary fields a{x) and tt{x) are defined on the dual 
lattice sites. The symbol < x,x > represents the 16 dual sites x adjacent to the direct 
lattice site x, and e{x) is the alternating phase (— 1)^0+^1+^2+^3 _ ggg [|iq| ^qj- more detailed 
discussion of the auxiliary field formulation used here. P^^ denotes the conventional Wilson 
plaquettes: 

The separate gauge couplings (3(j = Pt = Q/dr ^^^^ the extra parameter ■jp are 

introduced in the usual way to maintain Euclidean invariance in the continuum limit, when 
studying the model on anisotropic lattices at nonzero temperature. 

This paper studies the asymptotic scaling of the xQCD model and its perturbative 
thermodynamics. Section |I| studies how the four-fermion coupling affects the approach 
to the continuum limit and the lattice scale parameter Al. Section |T| deals with the 
perturbative thermodynamics of the model and presents the one-loop gluonic corrections 
to the energy, entropy, and pressure. Throughout the paper, the four-fermion interactions 
are dealt with via a mean field (or, equivalently, Nf 00) approximation, the details of 
which are given in Sec. Treating the model in this way, the effects of the four-fermion 
interactions on the asymptotic scaling behavior and the perturbative thermodynamics are 
governed solely by the value of the dynamical quark mass = (a)^ -|- (tt)^. Therefore, the 
one-loop gluonic corrections to the model's thermodynamic observables reduce to those of 
standard lattice QCD with a nonzero bare quark mass. Since such calculations have only 
appeared in the literature before at zero quark mass |1TT1-[13|, they are extended to nonzero 
masses in Sec. |V]. 



II. xQCD SCALE PARAMETER 



Consider the xQCD model on an isotropic lattice {ar = a), with /?o- = Pr = G/g^ and 
'jF = 1- In the asymptotic scaling region of lattice QCD the gauge coupling g is related to 
the lattice spacing a by 



aAL = Wog')-''^''?^ 



exp 



2Pog' 



l+0{g' 



(5) 



which defines the conventional lattice QCD scale parameter Al. In the continuum limit the 
scale parameter is the only dimensional parameter. Hence, all physical quantities measured 
on the lattice (such as particle masses, string tension, etc.) must be expressible in terms of 
Al by simple dimensional analysis. 

The numerical value of Al depends upon the specifics of the lattice regularization, and 
so is not universal. The ratio of the lattice and continuum scale parameters, in any two 
specific regularization schemes, is given by 



3 



— ^ = Maexp 
A. 



1 / 1 



(6) 



where M is the continuum renormahzation mass scale of interest (momentum subtraction 
point or Pauh-Villars mass, etc.). This ratio is most easily calculated by a one-loop back- 
ground field calculation 
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where the gluonic action is expanded around a slowly varying 
classical background configuration. 

In xQCD one has the additional four-fermion coupling G, which has physical dimensions 
of (length)^. In the continuum limit such an interaction is perturbatively irrelevant. There- 
fore, the asymptotic scaling formula (§) must remain valid in some window of a — * for any 
given value of G.^ The numerical value of the scale parameter will change depending on the 
particular value of G used. However, this dependence must cancel out when any physical 
quantities are expressed in terms of the continuum QCD scale parameter. This will need to 
be verified in future Monte Carlo simulations. 

The effects of the explicit four-fermion interaction in xQCD are most easily observed by 
a mean field calculation, or a large fiavor expansion, where the induced dynamical quark 
mass is related to the mean values of the auxiliary fields (see Sec. The calculation of the 
lattice scale parameter in lattice perturbation theory is then only affected by the presence 
of this dynamical mass in the fermionic contributions to the gluon self-energy (vacuum 
polarization). 

The relevant diagrams are: 



4 4 Jq/2 



Tilj^Spiq - k/2)-i^SF{q + k/2)] cos(g^a^) co?,{qua^ 



(7) 



Nfl r 

~4 Jq/2 '^^t^'^'^^^^)]^"/^ sin(g^a^)5^^. 



(8) 



(The notation is explained in the appendix.) Following |[T^|T^ and collecting together the 
gauge, ghost, and fermion contributions, the ratio of the xQCD scale parameter to the 
scale parameter Amin in the minimal subtraction scheme can be written 



A. 



A 



exp 



MIN 



/5oU8 2 ' 



(9) 



with 



^ Unless, of course, their exists a second order phase transition at some value of G = Gc, where 
the four-fermion interaction becomes both renormalizable and nontrivial. In this case, Eq. (|5|) may 
only hold when G is taken smaller than Gc- 
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J = i(ln47r - 7Euier) = 0.9769042, 
P = 0.0849780, 

f d'^q cos^(gi) cos^(g2) - I cos(2gi) cos(2g2) 



P. 



f (2vr) 



2 d'^q 
~ 3 J-oo (27r)^ 



(g^ + m?a?Y (g^ + 1)^ 



(10) 
(11) 



(12) 



Here, = (o")^+(7r)^ is a function of the two bare couplings {G, g}. In simulations of xQCD 
the value of the four-fermion coupling is to be chosen very small, such that the observed 
values of (a) and (tt) are principally due to nonperturbative effects of the gluons. Therefore, 
they cannot be reliably calculated in perturbation theory and must be determined by Monte 
Carlo simulation. For this reason, it is best to keep parameterized by the dynamical 

quark mass ma = a^J (cx)^ + (77)^. 

Values of Aum/A^ are plotted as a function of ma in Figs. |I| and || for Nj = 2 and 
Nf = 4, respectively. At ma = (G = 0) I find 



A 



MIN 



A 



MIN 



A. 



At 



10.846 Nf = 
16.518 Nf = 2 
28.779 Nf = A 



(13) 



in agreement with |jT6[. As ma increases the ratio increases monotonically. Assuming rea- 
sonable values such as (ipip) < 1 and 7 > 10, ma will be 0.1 or less. At ma = 0.1 the 
increase in Aum/A^ is about 0.15% (0.3%) for 2 (4) fermion flavors. Thus, it is a relatively 
small effect which will be difficult to extract from simulation data. 



III. xQCD THERMODYNAMICS 



Working in the canonical ensemble, the Helmholtz free energy F = —TlnZ, where Z = 
Trexp {—S/T) is the lattice regularized partition function. The temperature T = l/aT-N^-, 
where Nr (A^o-) is the number of lattice sites in the temporal (spatial) direction(s). Following 
the usual convention define ^ = a/ar- The energy density and pressure are then given by 



P 



T- 



1 dlnZ 
V 

d\nZ 



T 



dV 



1 T 

'S{aN^) 



dS 



rp4 



da J ' \iN, 
dS\ _ / Nr ^3 

'9^/ ~ ~Y\Jn^ 



\ 



Using the thermodynamic relations F 
can be expressed as 



U — TS and p 



,dS\ 

^d^/^Vda/\ 
-{dF I dV)T the entropy density 



dS\ 
I. — ) 



(14) 
(15) 



e + p 



4 e 
3T 



dS\ 
da r 



T 3T 3 V^A^^. 
The deviation from ideal gas behavior is typically characterized by 



(16) 
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Before proceeding any further, note that Gaussian integration over the auxihary fields 
appearing in Eq. (|l]) will yield a factor of {STc^/Nf'ya^)^'^^^ in addition to the desired four- 
fermion interactions. Since the thermodynamic observables are to be computed by differen- 
tiating the partition function with respect to a and ^, it is best to remove this factor from 
Z. This can be accomplished by adding to the action the extra term 

l^J2^^^rf]{x)r]{x), (18) 

where f] and t] free Grassmannian fields. Integration over f] and rj will then cancel the 
cancel the corresponding contribution from the auxiliary fields. Of course, if one is only 
concerned with measuring quantities relative to the zero temperature vacuum this becomes 
unnecessary. 

Adding (|T8|) to the xQCD action ([^) it follows immediately that 



+ «^^(XA<.X>+ (-7 + ?|)^((.^ + .=> + (.-„))]. (19) 



All the averages on the right hand side are per space-time volume. (Pss) {{Pst)) stands for 
the average value of the space-space (space-time) plaquettes. Note that {xQx) = ~3 

= -1- 

In order to calculate the regularization dependence of the bare parameters appearing 
above, a mean field treatment of the auxiliary fields together with perturbative QCD cor- 
rections will be used. The details are presented in Sees. |V| and 0. In this approach the 
calculations of d[3„/d^, dPr/d^, and d'jp/d^ become the same as in ordinary lattice QCD 
with (cr)^ + (tt)^ substituting for the square of the bare quark mass. The quantity d-y/d^ is 
determined by absorbing the ^ dependence of the one-loop gluonic corrections to the aipil) 
vertex into the value of G = I/7. This insures that the dynamical fermion mass, given by 
the pole in the quark propagator, is invariant with respect to ^. 

In Sec. |IV| it is shown that 



The quantity {dCm/d$,)^=i is plotted in Fig. as a function of ma = a^J (a)^ + (Tr)^. It is 
the only quantity of interest which varies significantly with ma. [{dCm/dC,)^=i varies about 
2% from ma = to 0.1, while dCa-/d$,, dCr/d^ and dCp/d^ all vary 0.2% or less over the 
same range.] Hence, the values of dp^/d^, dPr/d^ and d'-yp/d^ calculated with massless 
quarks can safely be substituted. 

Differentiating the action with respect to the spatial lattice spacing yields 



7 



/ dS_\ 
\ da I 



da 



da 



d\naj 4 



da 4 



+ 47 + a 



97 



da) 8 



AT 



(21) 



Using the identity ^(x(o"+^£^vr)x) = 2(1— 7^(cr^+7r^)) the last two terms can be combined. 
To lowest order in g, the values of dp^/da, d[3r/da and d'jp/da are the same as in ordinary 
lattice QCD, since they are controlled by the lowest universal coefficient of the beta function 
(see Sec. |Vp. There is no bare quark mass in the action. The expected scaling of such a 
mass is taken over here by the scaling of the four-fermion coupling. In Sec. |I^ it is shown 
that 



N 



dr 

7V daj^=i 



(22) 



Measuring the energy density and pressure with respect to the zero temperature vac- 
uum, by subtracting measurements on a symmetric lattice (denoted by (■ ■ ■)sym), the final 
expressions at ,^ = 1 are 



t^nU — H 



6 
9 

+ [l + 9 



9 



>dC^ 



d^ 

^dCF\Nf 



(Pss) - (Pss) 



sym 



di J A 
dCr, 



sym 



+ ( 1+^7^2:1-1^)7^ 



d^ 



sym 



T 

^ 3 3 



{Pss) ~ {Pss)sym + {Pst) ~ {Pst) 



sym 



27r2 



Nf 



7- 



sym 



(Pst) - (Pst) 



sym 



(23) 



(24) 



The (absolute) entropy density at ^ = 1 is 



4T3 



-N 



,lfdCr dC^\ 



(Pss) 



St 



2dCF\Nf 



2v' + ^2vr2 



9 1 9 dC^ 



di 



7^(a2 + 7r2)-l 



(25) 



In writing the above expressions I have used the identity {x^x) = — (TrAQ ^). The relations 
(|66|) and (^) from Sec. ^ have also been used to simplify the expression for the entropy 



density. For ma = ay (o")^ + (Tr)^ less than 0.1 the values of C^, Cr, and Cp are effectively 
the same as in ordinary QCD: 



V d^ 



Nf 

0.20161 - -^0.00062, 



(26) 
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0.13195 - ^0.00782, (27) 

0.2132. (28) 

They are discussed in full detail in Sec. |V]. Values of {dCmdC)i=i are tabulated in Table | 
for values of ma < 0.2. 

Note that 7^(cr^ + vr^) approaches the Gaussian free field value of 1 as 7 cxd. In 
which case, Eqs. (|23D - ([25| ) yield the correct expressions for massless QCD. Note also that 
in this model (Trgg-^) = 3 and ^(Tr(a + ieTx)Q-^) = 2(7^(^2 + tt^) - 1). Therefore, 

{TTMoQ-')sym = I - wji^^i^' + ^')smn " !)• (29) 

As of this writing the only numerical study of xQCD is that of Ref. [0, with Nf = 2 
massless flavors of quarks on a 8^ x 4 lattice at 7 = 10. In that study T = measurements 
(on a 8^ X 24 lattice) were only made at G/g"^ = 5.4. At that value of the gauge coupling 
they found {TiMoQ-^) - | = -0.0200(3) and -f^{a^ + vr^) = 1.0198(2) fl^, in good 
agreement with Eq. ( p9D and the vanishing of the entropy at zero temperature. Further zero 
temperature measurements would be necessary to extract the energy density and pressure. 

Even in the absence of such zero temperature measurements the 8^ x 4 lattice data from 
can be used with Eq. (|^) to extract the entropy density. The results for the gluonic and 
fermionic entropy densities are plotted separately in Figs. |^ and ^, respectively. The upper 
lines are the tree level results and the lower lines are the one-loop corrected results. The 
figure shows a clear first order transition consistent with the measurements of the Wilson 
line and the chiral order parameter in the same study. Theoretically |jl]], one expects a 
second order transitions in the continuum theory. The first order transition here is most 
likely due to the small size of the lattice used. Similar results are found in conventional 
lattice QCD with staggered quarks on lattices of that size. Clearly, larger lattices will be 
necessary to determine the nature of the chiral phase transition and the equation of state of 
the quark-gluon plasma in the continuum limit. 



dCr 
dCp 



5=1 



IV. MEAN FIELD APPROXIMATION 



In this section, the mean field treatment of the four-fermion interactions is discussed 
in greater detail and then used to explicitly calculate the regularization dependence of the 
four-fermion coupling G = I/7 to lowest order in the gauge coupling g. For these purposes 
it is best to work with the model where the aipip vertex and the dynamical fermion mass 
are proportional to \fG. This is accomplished by rescaling the auxiliary fields such that the 
action is 

x '- i<j j a=l ^ x,y j=l 

+ T.nx)^ E Hx)+ze{x)7r{x)]x''{x)] + ^j:i^\x) + n\x)]. (30) 
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FIG. 3. Gluonic entropy density using the data from 0, for two massless flavors on a 8^ x 4 
lattice at 7 = 10. The upper line is for the tree level results and the lower line the one-loop 
corrected results. 
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FIG. 4. Fermionic entropy density using the data from Q, for two massless flavors on a 8^ x 4 
lattice at 7 = 10. The upper line is for the tree level results and the lower line the one-loop 
corrected results. 
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+ 
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FIG. 5. Contributions to (a) up through 0{g^). 



')) it is easy to verify that 



Using identities such as \fG^{x{(y + isTT)x) = 2(1 — ^(cr^ + tt 
the final expressions of the previous section, ( |23D - (^), will remain unchanged. 

In a mean field (or a large Nf) approximation of the four-fermion interactions the aux- 
iliary fields take on the values 



(vr) 



GTt[S'p{x - 



x)], 



X 



x)], 



where Sp is the full quark propagator (including the gluonic contributions): 
Sy^ip) = iXl7Asin(pAaA)/aA + {a) + ii-,T?,{'K) + S(p) + 0{a). 



(31) 
(32) 
(33) 
(34) 



(35) 



Here S(p) is the quark self-energy of Sec. [VB| evaluated at rn? = (cr)^ + (vr)^. The lowest 
order gluonic contributions to (a) are shown in Fig. |^. 

Now the idea of the xQCD model is to choose the value of G to be very small, so that 
the breaking of the chiral symmetry is due to the long distance behavior of the gauge fields 
(as opposed to the spontaneous symmetry breaking in Gross-Neveu fl^ and Nambu-Jona- 
Lasinio type models). Since this is inherently a non-perturbative effect it will not be 
found in perturbation theory at any finite order of the gauge coupling. Therefore, we will 
have to rely on the values of (cr) and (vr) being determined by Monte Carlo simulation. 

For simplicity of the following discussion choose (tt) = 0. A non-vanishing value of (a) 
then signals dynamical chiral symmetry breaking. Following the one-loop calculation of the 
quark propagator in Sec. |V B| , the dynamically generated quark mass (determined by the 
pole in Sp{p)) is then 



y/G{a) 



l+g'-m) + 0{g' 



(36) 



Hence, the regularization dependence of this physical mass can be absorbed into the bare 
four-fermion coupling G by taking 



= VG{a) 



l+g'C„M) + 0{g' 



(37) 



with CmiO evaluated at m = vG'(cr). 
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FIG. 6. Contributions to the aipil) vertex up through 0{g^ 



The ^ dependence of G can also be determined by calculating the one-loop gluonic 
contributions to the aipip vertex shown in Fig. |^. It suffices to consider external momenta 
less than 1/a. It is then easy to see that the one loop gluonic correction to F^.^^ is identical 
to the scalar part of the one loop fermion self energy calculated in Sec. |V lj|: 



Nf 



2N 



Vg 



/ ^f^^Fip + q)SFip + k + q)-ff,D{q) 

X cos((p/, + q^/2)a^) cos{{p^ + k^ + g^/2)a^) 
+ terms higher 0(pa, ka) 



4-|Ai(i: 



3 J, 2Ai(0(A2(0 + G'(a)2)/a4 



N 



f 



1+/-S2(0 



(38) 
(39) 

(40) 



where S2 is the same as for ordinary QCD [Eq. (^)] with m = \fG{a). Defining F^^^(^) 
Z^\OK%^ the e dependence is absorbed by taking = ^fG[l+g%^iAi)-Mi))] 

a/G[1 + g'^Cm{i% in agreement with Eq. (jS^). 
Differentiating this with respect to ^ yields 



5=1 



(41) 



The scaling of G with the spatial lattice spacing (at fixed S,) is then given by G{a) = 
G[l + gmj:i,a - ^2)div], where i^i,^)Div = -l/iSir^) ln(a) and (Sa)^/!/ = -l/i^n^) ln(a) 



(42) 



are the divergent parts of S(p) as a ^ |2^. This yields 



G~d^ 



V. LATTICE QCD THERMODYNAMICS WITH NONZERO QUARK MASSES 

This section briefly reviews the perturbative thermodynamics of lattice QCD, including 
the contributions due to nonzero bare quark masses. The necessary one-loop calculations of 



ri|^13[ are then extended to nonzero quark mass and the results discussed. 

On an anisotropic lattice the action for QCD with Nf (degenerate) flavors of staggered 
fermions is 
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Nf/A 

X '- i<j j a=l x,y 



(43) 



with 



(44) 



The separate gauge couphngs Pa = Q/gl, Pt = Q/dry &s well as the extra parameter jp, 
are necessary to maintain Euclidean invariance (i. e. regularization independence) in the 
continuum limit. The functional dependence of these parameters on ^ = a/a,- is fixed by 
requiring the theory to be independent of ^ in the continuum limit. 

As shown in Sec. |T|, to calculate thermodynamic quantities one needs 



drriQa N 



and 



dS_\ 
da / 



+ a 



4 
da 



4 



dZ,,\N 



J 4 



(45) 



da 



da 4 



3 + 



dZ^ 
9 In a 



(46) 



The averages appearing above can be evaluated by Monte Carlo simulation. (xQx) = 
—Tt{QQ~"^) = —3, for 3 colors of quarks. Hence, when subtracting T = measurements 
on symmetric lattices the last terms appearing in Eqs. 



and (J6) will not contribute. 
In general, the various derivatives of the coefficients may also need to be calculated non- 
perturbatively As a first step, however, they can be estimated by lattice perturbation 
theory. 

Considering first the ^ dependence: 



PaiO 



6 



9^ 



MO 



i + CaiOg' + oig' 
i + CAOg' + oig') 

mo\l + CUOg' + 0{g' 



(47) 
(48) 
(49) 
(50) 



[Here g = g^^riC, = 1) and rriQ = mo{C, = 1).] The coefficients C^, Cr, Cp and Cm all vanish 
at ^ = 1 where most Monte Carlo simulations are performed. However, their derivatives 
with respect to ^ are nonzero at ^ = 1. They are calculated in Sees. |V A| and |V Ij . 

The implicit scaling of the gauge coupling g and the bare quark mass itiq with the spatial 
lattice spacing is governed by the renormalization group equations: 
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-a^ = P{g) = -Pog'-Pig' + 0{g') 







11 / 2Nf 
1 ^ 



167r2 
1 

7o 



33 y 167r2 

2 , ml 



102- ( 10 + -jATy 



27r 



2 ■ 



(51) 
(52) 
(53) 
(54) 



is universal (regularization independent) up through ©((yf^), and 7(5') is universal 
through 0{g^). At ^ = 1 



a— — = a 



da 



da 



d{6/g^ 
da 



-12/3o + 0(/), 



and 



oa 



(55) 



(56) 



The final perturbative expressions for the energy, pressure, and entropy (at ^ = 1) are: 



6 



r 1 9 
\_g2 



>9a 



e = T'N^l^il-g'^]{P,s) 



6 



+ moag 



^^^-^^(Trg-^) + 3^(-l ' ""'^ 



4 



a^ 



+ moa(l + /7o)^(Trg-i) + 3^ (3 + |^ 
4 4 V a ma 



(57) 



(58) 



4T^^4 r 6 
3 ^ir 



6 
9 



dCr' 



^dCF\Nf , . , ^ 1, 3Nf f . (9Z^ 
+ m„a[ 1 + 9^0 + 9'4^')${Tr(3-') 



a In a 



9^ / 16 

^ , ^v-^a aa 



{Pss) - (Pst) 



2dCF\Nf 



{TrMoQ-') - ^ + ^moa(TrQ-i) 



(59) 



(60) 



The mo ~ relations given by Eqs. (|66D , (^) and ( ^31) have been used to simplify the 
expression for the entropy density. In most cases one is interested in measuring the energy 
density and pressure relative to the T = vacuum, by subtracting measurements on sym- 
metric lattices. In which case, the last terms in Eqs. (0) and ( pBD will drop out. Using 
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{TtQQ~^) = 3 it follows that (TrA4oQ~^)sj/m = | — \moa{TTQ~^) sym- Hence, it is easy to 
verify that the expression for the entropy vanishes on an isotropic, symmetric lattice. 

The one loop perturbative calculations of the various coefficients appearing above, and 
in particular their dependence on the bare quark mass mo, are presented in the following 
two sections. 



A. Vacuum Polarization 

Define C^iO = C^iO + (0 and a(0 = C^iO + Cf (0 where the superscripts G 
and F refer to the gauge and fermionic contributions to the QCD vacuum polarization, 
respectively. The ^ dependence is most easily calculated using a background field method 
i. e. by expanding the gauge action around a classical background gauge configuration 
and requiring the effective action to be independent of ^. 



The pure gauge contributions have been calculated by Karsch [|rT|. They yield 
(9C^/a0c=i = 0.20161 and {dC^ /d^)^=i = -0.13195. The fermionic contributions are 
given solely by the one-loop vacuum polarization graphs of Eqs. (|^) and (^). These have 
been evaluated previously by Trinchero for massless quarks |12|. Following his notation. 



with 



Uij = {-hk.jh + + klU), (61) 

n^o = {-Kko + 5^ok^)Ir, (62) 



_ _iV/ r cos^(gia) cos^(g2a) - \ cos(2gia) cos(2g2a) , . 

_ ^iV/ r cos^(gia) cos^(goa^) - | cos(2gia) cos(2goa^) 

-'tUJ ~ o / FA fC\ I ^2^212/^4 ' / 



2 Jg/2 [A^iO + mWV/a^ 

and are related to these integrals by 

= 2[/.(0 - /.(I)], = 2[IA0 - m]- (65) 

For mo = I find {80^/00^^=1 = -^0.000622 and ((9Cf /90€=i = -^0.007822, in 
agreement with []T2[. The behaviors of ((9Cj/9^)g=i and {00^/8^^=1 (for Nj = 2) as a 



function of moct are shown in Figs. ^ and §, respectively. As moa increases from zero to 0.1, 
8CJ8i = 8C^/8i + 8C^/8i decreases 0.1% and 8Cr/8i = 8C^/8i + 8C^ /8i increases 
0.05% (per Nf/2 fiavors). Therefore, in most cases when moa is small these deviations can 
be ignored and the massless quark values used. 

The invariance of the string tension [|11[] as well as the vanishing of the gluonic entropy 
at zero temperature requires 

8^,\ ^fd^\ l 8{Q/g^) _ ^gg^ 



(9^ /5=i \ 8S, J 2 9 In a 

To lowest order in g this reads (5Co-/9^)g=i + {8Cr/80^=i = Pq. This relation no longer 
holds at nonzero moct. However, given the relatively small variations in {8C^/8C,)^=i and 
{8Cr/ 80^=1 it is still effectively satisfied for small bare masses. 
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-0.0006 



-0.0009 



-0.0012 



-0.0015 r 1 

0.00 0.04 0.08 0.12 0.16 0.20 

ma 



FIG. 7. {dC^ /d^)^=i as a function of the mass (for two quark flavors). 



-0.0072 



-0.0074 



c: -0. 



0076 



-0.0078 



-0.0080 ^ 1 

0.00 0.04 0.08 0.12 0.16 0.20 



ma 



FIG. 8. {dCl^ /d^)^=i as a function of the mass (for two quark flavors). 
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B. Quark Propagator 



jriO can be found by calculating the one loop corrections to the fermion propagator 
(i.e. the fermionic self-energy) and then demanding rotational invariance of the propagator 



in the continuum limit. This has also been carried out by Karsch for massless fermions [|13 
and will be extended here to the more general case. 
Consider 



To lowest order in g, — is given by the sum of the two terms: 



4CU 



,N^-1 



9'^ I ^f^SpiP + qh,iD{q) cos2(p^a^ + g;,a^/2) 



(67) 



(68) 



N'^ — II r 
9^ 2iV 2^"^^^^^^"'-^^"^'' J 



(69) 



Expanding the self energy in powers of py, 



nv) = 9' 



2N 



+ terms of higher 0(pa), (70) 



with 

Si, 
and 



a'^sin (qia) 



,,2Ai(0(A2(0+mga2) 
a"^ sin^(goaT) 



1 3 — cos(gia) — iAi(l) 
2+ A^i^) 
1 2 3 - cos(goa^) - |Ai(l) 
2+^ Al(0 



4-|Ai(i: 



4Ai(0 



, 2Ai(0(A2(0 + rriga2)/a4- 



Hence, 



+ mo 



1 + 9' 



1 + 9' 



2N 
2N 



Si,.(0 



+ iloPo 



ip/i + 9' 



2N 



Sl,r(0 



Demanding rotational invariance of T^"^^ in the continuum limit yields 

M0=a^ + 9'CF{0 + Oig% 



(71) 
(72) 

(73) 



(74) 



(75) 
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^ F -0.2140 - 



-0.2145 



-0.2150 

0.00 0.04 0.08 0.12 0.16 0.20 

ma 



FIG. 9. {dCp / d£,)^=i as a function of the mass. 



with 



_ 1 

Ci^(0 = ^^(Si,.(0-Si,.(0). 



(76) 



For mo = 01 find {dCp/d^)^=i = —0.2132, in agreement with |[T^. For nonzero quark 
masses {dCF/d^)^=i is shown in Fig. It decreases less than 0.2% as rriQa varies from zero 
to 0.1. Thus, in most cases the mo = value can be safely used. 

The remaining ^ dependence of F*^^) can be absorbed by a wavefunction rescaling: 



symi 



2N 



and a functional bare mass dependence 



Ar2 



2N 



l + g'Cm{0 + O{g' 



Si,.(0-S2(e) ■ 



(77) 
(78) 

(79) 
(80) 



At mo = ((9Si_a-/9^)^=i = —0.0368 and decreases just 0.2% as mo increases to 0.1. The 
divergent part of Si o- at ^ = 1 is [Sio-(l)]D/y = — l/(87r^) ln(a) [^O]- Thus, dZ^/dlna = 

Turning to the mass. 
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FIG. 10. {dCm/d£,)^=i as a function of the mass. 



/dmo\ 



fnog 



,fdC„ 



(81) 



{dCm/d^)^=i = —0.066 at niQa = 0. Figure |T0| shows the behavior of {dCm/d^)^=i as a 
function of rriQa. Its value decreases about 2% as rriQa increases from zero to 0.1. This is 
significantly more variation than seen in dCp/dC,, dC^dS, and dCr/d^, and is due mainly to 
the variation in 8122/ d^- Values of {dCm/ 9^)^=1 are also given in Table | at selected values 

of TTlQa. 

Requiring the entropy to vanish on an isotropic, symmetric lattice yields 



9 



+ 4 



and 



-9 



+ 9'^ 



dC„ 



5=1 



9 In a 







+ 



din mo 
9 In a 



0, 



(82) 



(83) 



which are satisfied by the numerical values given for moa ~ 0. Like the corresponding sum 
rule for the gauge couplings (^) they are slightly violated at larger moa. 
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APPENDIX 

On an anisotropic lattice with ^ = a/a-r the free {g = 0) lattice propagators 

m = ^ ^ 

Sf{p) = 



1 -i Ea 7a sin(pA)/aA + m 



where 



1=1 

3 

1=1 

Integration limits have been denoted by 

f _ fa d'^q r _ d!^q 
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ma 


/^/ 


ma 




0.000 


-0.0660 


0.100 


-0.0674 


0.005 


-0.0660 


0.105 


-0.0676 


0.010 


-0.0660 


0.110 


-0.0677 


0.015 


-0.0660 


0.115 


-0.0678 


0.020 


-0.0661 


0.120 


-0 0680 


0.025 


-0.0661 


0.125 


-0.0681 


0.030 


-0.0662 


0.130 


-0.0682 


0.035 


-0.0662 


0.135 


-0.0684 


0.040 


-0.0663 


0.140 


-0.0685 


0.045 


-0.0664 


0.145 


-0.0687 


0.050 


-0.0664 


0.150 


-0.0688 


0.055 


-0.0665 


0.155 


-0.0690 


0.060 


-0.0666 


0.160 


-0.0691 


0.065 


-0.0667 


0.165 


-0.0693 


0.070 


-0.0668 


0.170 


-0.0694 


0.075 


-0.0669 


0.175 


-0.0696 


0.080 


-0.0670 


0.180 


-0.0697 


0.085 


-0.0671 


0.185 


-0.0699 


0.090 


-0.0672 


0.190 


-0.0701 


0.095 


-0.0673 


0.195 


-0.0702 



TABLE I. {dCm/d$,)^=i at selected values of ma. 
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